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Wall  pressure  fluctuations  generated  within  rectangular  partial  enclosures  (RPEs)  have  been 
studied  experimentally  for  a  broad  range  of  geometrical  parameters.  The  geometry 
represents  a  simplified  version  of  a  new  generation  trapped  vortex  combustor  and  consisted 
of  a  rectangular  cavity  connected  to  a  neck  of  smaller  size.  Wall  pressure  fluctuations  have 
been  measured  through  wall  mounted  microphones  providing  single  and  multi-variate 
pressure  statistics  both  in  the  physical  space  and  in  the  Fourier  domain.  In  order  to  interpret 
the  pressure  signals,  aerodynamic  and  acoustic  investigations  have  been  carried  out  as  well 
for  several  cavity-neck  ratios.  The  analysis  of  the  acoustic  response  of  the  cavity  has  been 
conducted  both  numerically  and  experimentally  and  a  simple  theoretical  model  has  been 
proposed  to  predict  the  frequency  of  the  acoustic  resonances.  The  aerodynamic  study  has 
been  carried  out  through  PIV  measurements  that  provided  characterization  in  terms  of  the 
geometrical  parameters  of  both  the  large-scale  vortex  generated  within  the  cavity  and  the 
recirculation  zone  formed  upstream  of  the  neck.  The  use  of  the  POD  decomposition 
permitted  us  to  correlate  the  dynamics  of  the  recirculation  with  the  observed  pressure 
statistics.  The  aerodynamic  and  acoustic  investigations  allowed  us  to  interpret  exhaustively 
the  wall  pressure  cross-statistics  and  to  separate  contributions  induced  by  hydrodynamic 
and  purely  acoustic  pressure  fluctuations. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Rectangular  partial  enclosures  (hereafter  denoted  as  RPEs),  also  named  rectangular  covered  cavities,  are  encountered 
in  many  industrial  applications  including  heating,  ventilation  and  air  conditioning.  In  many  industrial  systems,  RPEs  are 
installed  to  control  heat  losses,  to  provide  convective  flow  partitioning  and  to  suppress  aerodynamic  noise  [1].  The  main 
drawback  in  the  use  of  RPEs  is  the  generation  of  pressure  oscillations  that  might  induce  structural  vibrations  and 
consequent  mechanical  failures  [2  .  Pressure  oscillations  generated  by  the  coupling  of  a  flow  with  acoustic  standing  waves 
are  commonly  called  self-sustained  or  self-excited  oscillations.  This  phenomenon  was  widely  investigated  in  the  literature  in 
T-joints,  safety  valves,  pipes  with  side  branches  and  open  cavities.  Typically,  the  interaction  of  a  shear  layer  with  the  walls 
has  been  identified  as  the  main  source  of  acoustic  waves  [3-6].  Several  observations  have  suggested  that  the  pressure 
oscillations  induced  by  aerodynamic  noise  into  RPEs  generally  fall  into  two  categories,  depending  on  whether  the  incoming 
jet  flow  is  free  or  bounded  on  one  side  [7-9],  the  latter  category  being  the  most  frequent  in  practice. 

Recently,  the  interest  for  RPEs  has  been  renewed  thanks  to  their  application  in  the  field  of  combustion  in  particular  in  the 
design  of  the  so-called  trapped  vortex  combustors  (TVCs)  10-12  .  In  this  combustor,  a  vortex  is  trapped  into  a  cavity,  this 
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configuration  being  beneficial  for  the  generation  of  a  stable  flame  and  an  efficient  combustion.  Despite  combustion  stability 
being  cited  as  a  relevant  characteristic  of  TVC,  Hsu  et  al.  [13]  identified  a  combustion  instability  in  which  the  flame  blow-out 
was  preceded  by  flame  intermittency  and  tonal  noise  emission.  More  recently  Agarwal  et  al.  [14  pointed  out  that  the 
acoustic  depth-mode  of  the  cavity  is  dominant  in  this  kind  of  combustors  and  it  influences  the  heat  release.  Gutmark  et  al. 
[15]  notice  that  TVC  is  susceptible  to  the  onset  of  combustion  oscillations  in  pressure  and  heat  release.  They  have  shown  that 
OH*  chemiluminescence  fluctuations  and  pressure  are  tonal  and  isofrequential,  for  specific  operative  conditions. 

Understanding  and  predicting  acoustic  instability  in  combustion  chamber  requires  the  knowledge  of  the  acoustic 
behavior  of  the  geometry.  Despite  the  aerodynamic  behavior  of  RPEs  being  studied  extensively  in  the  literature,  much  lesser 
studies  have  been  focalized  on  the  characterization  of  the  pressure  field,  and  of  the  statistics  of  the  wall  pressure 
fluctuations.  The  main  objective  of  the  research  presented  herein  is  to  provide  a  complete  description  of  the  pressure 
fluctuations  induced  in  a  RPE,  representing  a  simplified  version  of  a  real  TVC  in  non-reactive  conditions.  More  specifically, 
significant  geometrical  simplifications  have  been  applied  to  the  transverse  direction  to  make  the  resulting  flow  treatable  by 
planar  Particle  Image  Velocimetry.  The  resulting  geometry  enables  a  fundamental  study  in  the  flow  instability  mechanism 
affecting  a  spanwise-constrained  flow.  In  fact  a  more  important  three-dimensional  motion  is  expected  both  in  large  width- 
to-depth  ratio  and  in  axisymmetric  combustor  geometries. 

The  problem  was  studied  experimentally  and  for  several  geometrical  parameters.  Velocity  fields  were  determined 
through  the  Particle  Image  Velocimetry  (PIV)  technique  providing  an  overall  picture  of  the  aerodynamic  behavior.  The 
acoustic  response  of  the  cavity  was  analyzed  both  numerically  and  experimentally,  whereas  the  wall  pressure  fluctuations  in 
the  presence  of  the  mean  flow  were  measured  using  wall  mounted  microphone  arrays.  The  post-processing  data  have 
included  the  application  of  the  proper  orthogonal  decomposition  (POD)  to  the  PIV  fields  and  the  estimation  of  the  cross¬ 
statistics  (cross-correlations  and  cross-spectra)  of  the  wall  pressure  fluctuations. 

A  description  of  the  experimental  set-ups  investigated  is  given  in  Section  2,  while  the  preliminary  aerodynamic  and 
acoustic  characterization  of  the  RPE  are  presented  in  Section  3  and  4  respectively.  Section  5  includes  the  main  results  mainly 
consisting  of  the  wall  pressure  fluctuation  statistics,  whereas  conclusions  are  presented  in  Section  6. 


2.  Experimental  set-up 

2.1.  Test  case  and  flow  conditions 

Several  experiments  have  been  performed  for  obtaining  a  phenomenological  overview  of  the  self-excited  pressure 
oscillations  in  a  RPE.  In  Fig.  1  a  sketch  of  the  test  case  is  presented.  The  analyzed  geometry  represents  a  simplified  version  of 
a  so-called  third  generation  Trapped  Vortex  Combustors  that  are  typically  characterized  by  a  large  deep  cavity,  and  air  or 
fuel  injectors  positioned  on  one  side  of  it  [12,14].  In  order  to  modify  the  main  geometrical  parameters,  the  cavity  length,  L, 
was  kept  fixed  whereas  the  depth,  D,  was  varied  by  inserting  plates  on  the  cavity  bottom,  so  that  the  cavity  aspect  ratio  can 
be  defined  as  r  =  D/d ,  where  d  is  the  transversal  length  of  the  neck.  Although  in  the  present  study,  r  could  be  changed  from 
2.00  to  7.3,  particular  attention  has  been  given  to  three  geometrical  conditions:  r  =  3.8,  also  named  square  cavity ,  r  =  3.2 
named  rectangular  cavity  A  and  r  =  2.6,  named  rectangular  cavity  B  (see  Table  1).  As  illustrated  in  the  following,  the  cavity 
aspect  ratio  plays  a  fundamental  role  in  terms  of  flow-structure,  pressure  gradient,  acoustic  and  aeroacoustic  behavior.  The 
number  of  the  achievable  geometrical  configurations  is  72,  but  the  following  discussions  are  restricted  to  the  most 
interesting  cases.  In  Table  2  all  values  assumed  by  the  geometrical  parameters  are  summarized.  The  geometry  analyzed  is 
characterized  by  a  width  (W)-to-length  (L)  ratio  equal  to  0.33.  The  flow  in  the  narrow-span  rectangular  cavities,  W/L  <  1,  is 
two-dimensional.  This  assumption,  according  to  the  results  carried  out  by  16],  has  been  further  confirmed  by  several  RANS 
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Fig.  1.  Sketch  of  the  frontal  and  the  transverse  view  of  the  test  case  and  enlargement  of  the  inlet  section.  The  main  geometrical  parameters  are  also 
reported  for  clarity. 
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Table  1 

Sketch  of  three  RPEs  on  which  the  experimental  investigation  is  focused,  for  all  of  these  configurations  W=60  mm  and  d  =  50  mm. 


Configuration 

r 

D  (mm) 

Geometry 

Square  cavity 

3.8 

190 

_ 

Rectangular  cavity  A  3.2  160 


Rectangular  cavity  B  2.6  130 


Table  2 

Range  of  values  of  the  test  case  geometrical  parameters. 


r 

D  (mm) 

L  (mm) 

W  (mm) 

/  (mm) 

d  (mm) 

din  (mm) 

2-7.3 

100-220 

190 

60 

50 

30-50 

7 

and  URANS  simulations  that  have  been  carried  out  by  the  authors  and  that  are  not  reported  here  for  brevity.  As  rendered  in 
Fig.  1,  the  air  is  injected  into  the  test  case  through  the  three  circular  inlets  that  are  7  mm  in  diameter  (din).  The  jet  mass-flow 
rates  range  from  7  g/s  to  15  g/s,  for  a  corresponding  bulk  velocity  inlet  of  50  and  108  m/s.  Taking  into  account  the  bulk  inlet 
velocity,  Uin,  the  jet  Reynolds  number  follows 


Re  = 


Uind 


ln^in 
V  5 


for  the  mass-flow  values  under  investigation,  the  Reynolds  number  interval  is  2  x  104  <Re<5  x  104.  These  Re  numbers 
ensure  low-frequency  acoustic  wave  emissions  within  the  region  1.2  and  2.5  jet  diameters  [17,18  .  Such  perturbations  excite 
the  RPE  first  modes,  ensuring  the  achievement  of  the  self-excited  oscillations. 


2.2.  Air  supply  system 

The  facility  investigated  herein  has  been  designed  to  carry  out  experiments  at  atmospheric  pressure.  The  air  supply 
consists  of  two  lines:  the  first  one  connects  an  high  pressure  plenum  chamber  to  the  test  rig  through  a  series  of  ball  valves, 
pressure  regulators  and  flow  meters.  The  second  one,  denoted  as  bypass  line,  is  connected  to  the  PIV  seeding  system 
described  in  the  next  paragraph.  The  air  supply  system  is  depicted  in  Fig.  2  by  a  single-line  diagram.  The  flow-meter  adopted 
to  regulate  the  mass-flow  is  a  Brooks  Smart  Mass  Flow  5853 S.  The  flow-meter  is  connected  to  a  home-made  muffler- 
collector  very  similar  to  the  one  proposed  by  19].  It  is  designed  in  a  single-inlet/triple-outlet  configuration,  to  split  the  main 
mass-flow  in  three  equivalent  flows.  The  muffler-collector  is  filled  with  an  acoustic  absorbing  foam  in  order  to  damp  the 
narrow-band  tonal  noise  generated  in  the  air  supply  system  (see  Fig.  3).  Several  measurements  have  been  performed  to 
verify  the  homogeneity  of  the  mass-flow  at  the  outlets  of  the  collector  and  to  verify  the  effectiveness  of  the  acoustic 
damping  provided  by  the  muffler.  For  the  sake  of  conciseness,  these  results,  satisfactory  in  all  cases,  are  not  reported  herein. 

2.3.  Particle  Image  Velocimetry 

The  cavity  investigated  is  equipped  with  three  optical  accesses,  the  top,  bottom  and  side  walls,  suitable  for  the 
application  of  the  PIV  measurement  technique.  A  scheme  of  the  PIV  set-up  is  given  in  Fig.  4(a).  The  PIV  system  consists  of  a 
PCO  Sensicam  (double  shutter  camera,  1280  x  1024  resolution)  and  a  2  mm  thick  laser  sheet  created  by  a  double  cavity  Nd: 
YAG  laser  (Big  Sky  Laser/Quantel  200  mj  light  pulses).  Each  pair  of  acquired  images  has  been  processed  using  the  software 
PlVdef,  developed  by  researchers  of  the  Italian  Ship  Model  Basin  -  CNR  center  [20  .  The  post-processing  procedure  uses 
iterative  adaptive  correlations  and  windows  deformation  algorithms.  The  resolution  of  the  interrogation  area  has  been 
varied  from  64  x  64  pixels  to  16  x  16  pixels  with  50  percent  overlap.  After  each  iteration,  the  vector  map  has  been  filtered  to 
remove  spurious  vectors,  identified  by  median  tests  and  replaced  using  2D  interpolations.  Eventually,  a  peak  locking 
analysis,  based  on  the  displacement  probability  density  function  (PDF),  has  confirmed  the  proper  setting  of  the  final 
iteration.  The  time  delay  between  two  consecutive  laser  pulses  has  been  adjusted  from  50  ps  to  150  ps  to  provide  a  mean 
displacement  vector  of  about  5  pixels  and  a  velocity  uncertainty  of  1-2  percent.  This  accuracy  is  referred  to  all  flow 
conditions  investigated.  The  image  pair  acquisition  rate  was  5  Hz.  A  seeding  spray  of  DEHS  (DiEthylHexylSebacate) 
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Fig.  2.  Single-line  diagram  of  the  feedline  system. 


Fig.  4.  Three  dimensional  scheme  of  the  PIV  experimental  set-up  (a).  Sketch  of  the  three  different  FOV  investigated  (b). 


characterized  by  1  y im  diameter  particles  has  been  produced  by  a  PivTec  Laskin  nozzle.  The  seeding  was  mixed  with  the 
main  stream  air  and  subsequently  injected  into  the  muffler  as  shown  in  Fig.  2.  Essentially,  three  fields  of  view  (FOV)  have 
been  investigated:  the  first  one,  denoted  as  FOVt,  is  190  x  190  mm2  wide  and  covers  almost  the  whole  cavity;  the  second 
one,  named  FOV2,  is  90  x  90  mm2  wide  and  is  centered  on  the  neck;  the  last  one,  denoted  as  FOV3  and  60  x  60  mm2  wide,  is 
centered  on  the  top  wall  of  the  RPE.  A  scheme  of  the  three  FOV  is  reported  in  Fig.  4(b). 
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Fig.  5.  Sketch  of  the  static  pressure  taps  positioned  on  the  top  wall  of  the  RPE. 


Fig.  6.  Sketch  of  the  experimental  set-up  adopted  for  the  aeroacoustic  measurements. 

2.4.  Static  wall  pressure  measurements 

The  symmetry  plane  of  the  RPE  model  has  been  equipped  with  23  static  pressure  taps  0.9  mm  in  diameter,  placed  along 
the  top  wall.  The  array  begins  10  mm  downstream  of  the  RPE  backward  face,  and  the  streamwise  separation  between  two 
consecutive  taps  was  10  mm  (measured  from  center  to  center)  as  reported  in  Fig.  5.  An  Airflow  TSI  PVM620  series  pressure 
transducer  was  used  to  measure  a  differential  pressure  in  the  range  -3735  Pa  to  3735  Pa  and  a  rotary  scanivalve  was 
implied  to  sequence  the  pressure  ports. 


2.5.  Aeroacoustic  measurements 

Aeroacoustic  measurements  have  been  performed  for  several  mass  flow  inlet.  Five  microphones  flush  mounted  at  wall 
have  been  installed  as  illustrated  in  Fig.  6.  Microtech  Gefell  M360  microphones  have  been  used,  four  of  which  flush- 
mounted  on  the  RPE  and  one  installed  stand-alone,  100  mm  far  from  the  outlet.  The  microphone  outputs  were  acquired 
using  a  National  Instruments  A/D  board  at  20  kHz  for  60  s  using  an  anti-aliasing  cut-off  frequency  filter  at  5  kHz.  The  mass 
flow  inlet  was  varied  between  7  g/s  and  15  g/s,  corresponding  to  a  Re  based  on  the  inlet  diameter  ranging  between  about 
2  x  104  and  5  x  104.  The  aspect  ratio  was  varied  from  2.0  to  4.4  in  order  to  obtain  a  total  number  of  36  more  interesting 
geometries. 
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Fig.  7.  Sketch  of  the  experimental  set-up  adopted  for  the  modal  test. 


2.6.  Acoustic  measurements 

In  a  stand-alone  configuration,  in  still  air,  the  acoustic  properties  of  the  RPE  have  been  preliminary  tested  through  a 
modal  test  analysis  carried  out  experimentally  in  a  semi-anechoic  chamber  (see  [21]  for  more  details  on  the  anechoic 
chamber).  For  the  purpose  of  damping  the  vibration  of  the  system,  the  model  was  supported  by  a  tripod  positioned  on  a 
sand-box.  A  loudspeaker,  preliminarily  characterized,  was  located  150  mm  far  from  the  outlet  section,  with  the  scope  of 
forcing  the  acoustic  modes  of  the  test  case.  A  sketch  of  this  experimental  set-up  is  shown  in  Fig.  7.  The  modal  test  has  been 
performed  by  applying  a  swept-sine  spanning  from  100  Hz  to  500  Hz  and  using  the  same  measurement  chain  and 
acquisition  parameters  as  the  ones  adopted  for  the  described  aeroacoustic  measurements.  As  will  be  shown  below,  the 
spectral  response  is  given  in  terms  of  the  frequency  response  function  (FRFSPL)  that  is  computed  as  follows: 


FRFspl  =  (2) 

^ps 

where  LpS  is  the  so-called  reference  sound-pressure  spectrum  level  [22]  computed  from  external  reference  microphone 
located  outside  of  the  test  case,  inside  the  anechoic  chamber  and  in  close  proximity  of  the  loudspeaker.  Instead  Llps  is 
determined  inside  the  cavity  using  mic.  5.  Varying  D  from  100  mm  to  220  mm  and  d  from  50  mm  to  30  mm,  36  different 
cavity  aspect  ratios  have  been  tested  to  determine  accurately  the  dependence  of  the  first  mode  frequency  upon  f. 


3.  Fluid  dynamic  characterization 

The  aerodynamic  characterization  provides  an  overall  picture  of  the  main  flow  structures  forming  within  the  cavity. 
A  simplified  sketch  of  the  RPE  flow  physics  is  represented  in  Fig.  8.  For  all  the  aspect  ratios  analyzed  herein,  a  few  important 
flow  structures  have  been  identified:  a  free  jet  flow,  a  jet- wall  interaction  region,  a  large  scale  vortex,  a  recirculation  bubble 
formed  into  the  neck,  a  stagnation  point.  As  described  in  Section  2.3,  in  order  to  describe  with  the  proper  spatial  resolution 
all  these  flow  behaviors,  different  FOV  were  adopted  for  the  PIV  measurements.  Fig.  9(a)  and  (b)  reports  the  averaged 
velocity  field  for  a  fixed  mass  adduction  equal  to  15  g/s,  for  a  cavity  aspect  ratio  of  2.6  and  3.8  respectively.  It  can  be  observed 
that  the  fluid  dynamic  behavior  is  strongly  affected  by  r.  Indeed,  despite  the  r  variation  inducing  a  weak  modification  of  the 
vortex  displacement,  the  vortex  intensity  varies  significantly.  In  Fig.  9(c)  and  (d),  the  azimuthal,  V0,  and  radial,  Vr,  velocity 
profiles,  extracted  along  the  x  direction  using  the  vortex  center  as  a  reference,  are  reported.  It  should  be  stressed  that 
previous  studies  [23]  demonstrated  that  the  azimuthal  velocity  evolution  influences  the  mixing  and  the  main  vortex 
stability  into  TVCs.  It  is  noticeable  that  for  the  square  cavity  (7"=3.8),  the  Ve  distribution  is  symmetric  with  respect  to  the 
vortex  center  and  exhibits  a  linear  trend  in  analogy  with  the  inner  zone  of  a  Rankine  vortex.  This  result  suggests  that  the 
square  cavity  configuration  traps  the  vortex  in  the  most  efficient  way.  Indeed,  [11]  experimentally  investigated  a  geometry 
very  similar  to  the  present  one  but  in  reactive  condition,  and  concluded  that  the  square  cavity  configuration  is  the  most 
efficient  one  in  terms  of  combustion  performances.  Another  interesting  feature  consists  in  the  formation  of  a  recirculation 
bubble  in  the  neck.  This  effect  is  due  to  the  presence  of  an  adverse  pressure  gradient  and  it  is  observed  for  all  the  cases 
under  consideration.  The  length  scale  and  the  position  of  this  structure  are  affected  by  r,  since  this  parameter  influences  the 
adverse  pressure  gradient  intensity.  In  order  to  better  characterize  this  feature,  PIV  measurements  carried  out  in  the  FOV2 
have  been  analyzed.  Results  reported  in  Fig.  10(a)  and  (b),  corresponding  again  to  r=  2.6  and  3.8,  demonstrate  that  a 
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Fig.  9.  Averaged  velocity  fields  measured  at  Re  =  50  271  for  different  values  of  the  cavity  aspect  ratio:  r=  2.6  (a),  r= 3.8  (b).  Azimuthal  (o)  and  radial  (A) 
velocity  profiles  extracted  from  the  PIV  data  aty  =  73  mm  for  r= 2.6  (c)  and  aty  =  90  mm  for  r= 3.8  (d). 


recirculation  bubble  is  formed  in  the  neck  independent  of  the  aspect  ratio.  The  scale  and  the  position  of  the  recirculation 
region  are  strongly  affected  by  r,  its  size  being  larger  for  larger  r. 

In  order  to  better  characterize  the  influence  of  r  on  the  recirculation  bubble  intensity,  the  normalized  mean  static 
pressure,  measured  along  the  upper  wall  of  the  cavity  model,  has  been  analyzed.  The  results  achieved  for  r=  2.6  and  3.8  are 
reported  in  Fig.  11.  The  data  are  reported  in  term  of  Cp*  =  (pw  -pe) /JpUfn,  where  pw  is  the  time  averaged  wall  pressure,  pe  is 
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(a)  (b) 


Fig.  10.  Streamlines  of  the  averaged  velocity  field,  corresponding  to  FOV2  and  measured  at  Re  =  50  271  for  r=  2.6  (a)  and  r= 3.8  (b). 


Fig.  11.  Pressure  coefficient  distribution  along  the  top  wall  of  the  RPE:  (A)  square  cavity,  (o)  rectangular  cavity. 


a  reference  pressure,  measured  at  the  exit  of  the  model,  Uin  is  the  mean  flow  velocity  at  the  inlet  and  p  is  the  fluid  density. 
The  streamwise  distance  is  normalized  with  respect  to  the  top  wall  length,  i.e.  L+l.  It  can  be  observed  that  the  most  intense 
vortex  generated  for  r=3.8  induces  a  relevant  variation  of  pressure  in  the  cavity  and  thus  the  largest  adverse  pressure 
gradient.  In  summary,  a  significant  connection  among  the  large  scale  dominant  vortex  intensity,  the  recirculation  bubble 
properties  and  T,  is  evident.  The  results  to  be  presented  in  the  following  shall  demonstrate  that  the  cavity  aspect  ratio 
affects  significantly  the  wall  pressure  statistics  as  well  as  an  effect  of  its  influence  upon  the  topology  and  intensity  of  the 
large  scale  flow  structures. 

4.  Acoustic  results 

The  acoustic  modes  in  the  self-excited  condition  are  expected  to  influence  significantly  the  wall  pressure  fluctuations 
even  in  the  presence  of  a  mean  flow.  For  this  reason,  a  detailed  analysis  of  the  acoustic  response  of  the  cavity  has  been 
carried  out  with  the  scope  of  determining  quantitatively  the  typical  frequencies  of  the  acoustic  resonances.  A  theoretical 
model  has  been  developed  to  provide  a  prediction  of  the  frequency  of  the  first  acoustic  mode,  and  experiments  and  FEM 
numerical  simulations  have  been  performed  to  demonstrate  its  reliability.  Those  theoretical,  numerical  and  experimental 
results  are  presented  in  the  next  subsections. 

4.1.  Acoustic  numerical  simulation 

A  three-dimensional  finite  element  model  of  the  present  RPE,  using  the  commercial  software  COMSOL,  has  been 
developed  and  solved  in  the  frequency  domain.  The  pressure  distribution  is  provided  over  the  frequency  range  50-2000  Hz 
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Fig.  12.  Sketch  of  the  mesh  adopted  (a);  representation  of  the  boundary  conditions  and  the  probe  locations  (b). 


with  a  resolution  of  1  Hz.  All  sides  of  the  exterior  domain  surrounding  the  cavity  have  been  modelled  as  open  radiating  into 
the  free  space.  As  acoustic  forcing,  a  constant  pressure  has  been  imposed  at  the  vertical  wall  of  the  external  domain. 
Dissipative  effects  due  to  the  viscous  losses  of  the  air  moving  back  and  forth  in  the  neck  have  been  neglected.  A  mesh 
refinement  procedure  has  been  applied  in  order  to  achieve  solutions  independent  of  the  mesh  size  and  a  scheme  of  the 
mesh.  The  final  mesh  adopted  after  the  refinement  procedure  is  given  in  Fig.  12  (a).  This  configuration  corresponds  to  105 
tetrahedral  cells. 

The  boundary  conditions  assumed  for  the  simulation  and  the  probe  points  taken  into  account  to  compute  the  frequency 
response  function,  FRFSpl,  are  illustrated  in  Fig.  12(b).  FRFSpl  computed  from  the  numerical  simulations  reported  in  Fig.  13 
are  analyzed  to  identify  the  frequencies  of  the  modes  that  correspond  to  the  peaks  in  the  FRFSPL.  The  first  three  modes 
achieved  in  the  case  of  a  square  and  a  rectangular  A  cavity  are  reported  in  Figs.  14  and  15  respectively.  In  the  square  case, 
coalescence  of  the  frequencies  associated  to  the  second  longitudinal  and  transverse  modes  occurs.  As  a  consequence 
the  second  mode,  calculated  by  the  FEM  code,  appears  as  a  linear  combination  of  longitudinal  and  transversal  modes  in 
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Fig.  13.  Frequency  response  function  for  three  different  configurations:  r=  3.8  (— ),  r=  3.2  ( - ),  r= 2.6  ( . ). 


(a)  (b) 


Fig.  14.  First  three  modes  occurring  in  the  square  RPE.  Gray  scale  represents  the  sound  pressure  level  in  arbitrary  scale. 


Fig.  14(b).  This  combination  is  marked  as  2 LT.  The  values  of  the  modal  frequencies  obtained  numerically  are  summarized  in 
Table  3.  It  is  interesting  to  note  that  in  the  region  close  to  the  outlet  section,  a  development  of  the  first  mode  outside  of  the 
neck  can  be  clearly  observed  (see  Fig.  16).  Such  a  result  suggests  that  in  order  to  predict  correctly  the  mode  frequencies 
using  a  plane  wave  model,  an  appropriate  end  correction  has  to  be  applied.  This  approach  is  presented  in  the  next  section. 


42.  Plane  waves  model  and  experimental  results 

The  frequency  of  the  first  acoustic  mode  is  predicted  using  a  plane  wave  propagation  model.  The  plane  wave  assumption 
is  correct  only  for  low  frequency  oscillations,  a  reasonable  hypothesis  when  only  the  first  mode  is  taken  into  account.  As  a 
further  assumption,  the  medium  in  the  system  is  assumed  to  be  uniform  and  at  rest.  An  acoustic  lumped  model  can  be 
established  by  decomposing  the  RPE  into  two  regions  separated  by  a  section  discontinuity.  These  regions  are  commonly 
indicated  as  cavity  and  neck  (see  Fig.  17).  The  general  solution  of  the  homogeneous  wave  equation,  satisfying  the  hard-wall 
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Fig.  15.  First  three  modes  occurring  in  the  rectangular  RPE:  (a)  first  mode,  (b)  second  longitudinal  mode,  (c)  second  transversal  mode  and  (d)  third  mode. 
Gray  scale  represents  the  sound  pressure  level  in  arbitrary  scale. 


Table  3 

First  three  mode  frequencies  varying  r. 


r 

UT  (Hz) 

2 LT  (Hz) 

2 L  (Hz) 

2 T  (Hz) 

3 LT  (Hz) 

3.8 

185 

/ 

936 

936 

1273 

3.2 

202 

/ 

935 

1083 

1378 

2.6 

225 

942 

/ 

1321 

1534 

Fig.  16.  Details  of  the  exterior  acoustic  field  simulated  at  the  first  mode  frequency.  The  length  of  exterior  end  correction  5e  is  qualitatively  estimated  and 
reported  in  the  figure.  The  gray  scale  represents  the  sound  pressure  level  in  arbitrary  scale. 


and  open-end  boundary  conditions,  can  be  written  in  the  following  form: 

-L=  tan  (kL)  tan  (Id),  (3) 

where  k  is  the  wavenumber,  L  the  cavity  length  and  l  the  neck  length.  Eq.  (3)  was  found  by  24]  and  applied  by  [25]  to 
predict  the  frequency  of  the  first  acoustic  mode  of  a  so-called  double  duct  geometry.  A  similar  formulation  has  been 
suggested  by  26]  to  predict  the  first  acoustic  mode  of  backward  facing  step  geometries  in  combustors  (see  also  [27]). 
A  simplified  resonant  condition,  corresponding  to  kL  ^  1  and  kl<  1,  can  be  derived  from  Eq.  (3).  By  expanding  this  equation 
to  the  first  order,  the  following  conventional  expression  for  the  Helmholtz  resonator  (HR)  is  achieved: 

Y  =  k2Ll. 


(4) 
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Si=DW  S2=dW 


BC:  ui=0  BC:  p2=0 


Fig.  17.  Schematic  representation  of  the  lumped  model. 


Fig.  18.  Solutions  of  Eqs.  (4)  ( cubical  HR  — )  and  (3)  ( double  duct - )  obtained  with  Newton's  method. 


For  both  Eqs.  (3)  and  (4),  the  wavenumber  depends  on  the  neck  length,  the  cavity  length  and  the  cavity-neck  section  ratio,  i. 
e.  k  =f(T,U  0-  By  computing  L  and  l  for  the  present  test  model,  the  solution  of  Eq.  (3)  collapses  well  into  Eq.  (4).  In  Fig.  18  the 
departure  of  the  solution  of  Eq.  (3)  from  Eq.  (4)  is  reported  as  a  function  of  r.  It  is  worth  noting  that  despite  the  geometry 
under  investigation  being  quite  different  from  a  cubical  HR,  with  a  centered  and  small  orifice,  its  resonant  behavior  at  high  r 
can  be  represented  well  using  the  simple  HR  equation.  The  discrepancy  observed  at  small  r  is  related  to  the  fact  that  the 
present  simple  model  is  based  on  the  open-end  boundary  condition,  p2(L+t)  =  0.  As  was  shown  from  the  numerical  results 
in  Fig.  16,  this  assumption  is  not  representing  the  actual  physical  behavior  of  the  pressure  at  the  outlet  (see  also  [26  ).  As  a 
consequence,  in  order  to  take  into  account  the  dependence  upon  time  of  the  pressure  at  the  outlet,  an  appropriate  end- 
correction  has  to  be  introduced.  From  a  formal  viewpoint,  the  end-correction  is  commonly  denoted  as  follows: 

r  =  l+Se+5i  =  l+2SR,  (5) 

where  5e  and  Si  are  exterior  and  interior  end-corrections  respectively.  Therefore,  Eq.  (4)  can  be  modified  to  obtain  the 
following  expression: 

j  =  k2U,  (6) 

The  end-correction  term  commonly  adopted  in  the  literature  is  provided  by  the  following  expression: 


f 


=  l+5e  +  5i  =  1+28r  =  l+a 


where  y/Wd/x  is  the  hydraulic  radius  of  the  neck  [22  .  The  end-correction  for  axi-symmetric  HR,  expressed  in 
dimensionless  form,  depends  on  the  cavity-neck  hydraulic  diameters  ratio  and  it  ranges  between  two  values, 
0.785  <  2SR/d  <  0.848  [28].  It  has  been  checked  that  in  the  present  case,  these  literature  predictions  do  not  lead  to 
satisfactory  results  probably  due  to  the  peculiar  neck  shape  and  position  in  the  present  geometry.  Hence  a  novel  semi- 
empirical  estimation  of  SR  has  been  accomplished  on  the  basis  of  our  experimental  results. 

The  experimental  modal  test  performed  allows  us  to  calculate  FRFSpl  that  is  reported  in  Fig.  19.  The  frequency  and  the 
amplitude  of  the  first  mode  are  strongly  affected  by  r.  In  particular  the  amplitude  of  the  mode  increases  significantly  in  the 
square  cavity  case.  The  dependence  of  the  first  mode  frequencies  upon  r,  determined  by  means  of  the  experimental 
investigation,  carried  out  for  several  r,  is  reported  in  Fig.  20  in  a  log-log  scale.  In  this  figure  the  typical  Helmholtz  resonator 
power  law,  /  ~  T~1/2,  can  be  observed.  Hence,  despite  the  RPEs  investigated  being  quite  different  from  classical  HR,  it  can  be 
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Frequency  (Hz) 

Fig.  19.  Frequency  response  function  for  three  different  configurations:  r=  3.8  (— ),  r=  3.2  ( - )  and,  r= 2.6  ( . ). 


Fig.  20.  First  mode  frequency  referred  to  d=50  mm  (o),  d= 40  mm  (■),  d  =  30  mm  (o).  The  dashed  line  represents  a  power  law  decay  with  exponent  - 1/2. 


Frequency  ^  (Hz) 


Fig.  21.  Frequency  of  the  first  mode,  measured  experimentally,  versus  the  frequency  predicted  applying  the  empirical  end-correction  modelled  in  Eq.  (8). 
The  linear  fitting  of  the  data  is  superimposed  (— ). 
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concluded  that  the  first  acoustic  mode  is  exactly  an  HR  mode.  Whereas  the  weak  frequency  drifting,  which  occurs  as  a 
function  of  the  transverse  length  of  the  neck,  d,  is  due  to  the  development  of  the  first  mode  in  the  exterior  domain,  as  was 
shown  in  Fig.  16.  From  a  proper  fitting  of  the  experimental  data,  the  following  analytical  expression  of  a  is  derived: 


L  (2^1 01*)  1  \fwd’ 


where  q  is  the  intercept  of  the  linear  fitting  of  the  experimental  data  and  its  value  is  close  to  2.55  for  all  the  investigated 
cases.  A  comparison  between  the  experimental  frequencies  and  those  predicted  by  the  end-correction  proposed  in  Eq.  (8)  is 
reported  in  Fig.  21,  showing  a  very  good  agreement. 


5.  Wall  pressure  fluctuation  statistics 

In  the  presence  of  a  flow,  the  acoustic  pressure  previously  depicted  is  combined  with  the  complex  aeroacoustic  behavior 
generated  by  the  flow  structures  created  inside  the  cavity.  The  estimated  acoustic  mode  resonances  as  well  as  the  velocity 
fields  measured  with  the  PIV  technique  provide  a  valid  support  to  interpret  the  pressure  data  presented  in  both  the  physical 
and  the  Fourier  domain.  In  the  following  subsections,  single  and  multi-variate  statistics  are  also  presented.  The  achieved 
results  provide  an  overall,  reasonably  clear  picture  of  the  relationship  between  the  flow  structures,  the  acoustic  resonance 
and  the  wall  pressure  statistics,  parameterized  in  terms  of  r  and  Re.  We  refer  to  Fig.  6  for  the  microphones  numbering  and 
positions. 


5.1.  Single-point  statistics 

The  overall  amplitude  of  the  wall  pressure  fluctuations  has  been  characterized  through  the  computation  of  the  root  mean 
square  pressure  coefficient  Cprms,  defined  as  follows: 

Cprms=T%,  (9) 

ipUfn 

where  ov  represents  the  standard  deviation  of  the  pressure  signal.  As  shown  in  Fig.  22  the  largest  pressure  levels  are 
measured  in  the  cavity  (microphones  1,  2,  3  and  5),  whereas  an  important  reduction  occurs  in  the  neck  (microphone  4).  A 
quite  similar  behavior  is  obtained  for  each  Re  investigated.  The  relevant  reduction  of  Cprms  at  mic.  4  location  can  be  ascribed 
to  reflections  of  the  acoustic  waves  due  to  the  section  contraction  sensed  by  the  acoustic  energy  transmitted  from  the  cavity 
to  the  neck.  It  is  also  interesting  to  note  that,  except  for  r=  2.6  at  the  location  of  mic.  4,  the  Cprms  is  weakly  dependent  upon 
r.  In  order  to  better  clarify  this  behavior,  the  pressure  fluctuations  measured  by  mic.  4  are  analyzed  in  more  detail.  The  wall 
pressure  power  spectral  density  (PSD)  measured  by  mic.  4  is  reported  in  Fig.  23  for  a  fixed  Re  and  different  r.  Except  for  the 
case  at  the  lowest  7",  the  shape  of  the  spectra  is  similar  in  all  cases,  exhibiting  a  clear  peak  at  the  frequency  of  the  first 
acoustic  mode.  The  lowest  r  case  is  characterized  by  a  larger  energy  content  in  the  low  frequency  range  that  corresponds  to 
the  observed  increase  of  the  Cprms  magnitude.  This  behavior  can  be  ascribed  to  the  influence  of  the  adverse  pressure 
gradient  induced  in  the  neck  that  becomes  significant  in  the  lowest  r  case  (see  Fig.  11  for  0.8  <x/(L+t)  <  1).  This  aspect  is 
further  clarified  by  the  analysis  of  the  probability  density  functions  (PDFs)  of  the  pressure  fluctuations  measured  by  mic.  4 
that  are  reported  in  Fig.  24  for  different  aspect  ratios.  The  random  variable  is  represented  in  reduced  form  in  order  to  have 


Fig.  22.  Cprms  referred  to  all  microphones  and  obtained  for  Re  =  50  271  and  different  aspect  ratios:  r= 3.8  (□),  r= 3.2  (o)  and  r=  2.6  (A).  The  arrow  evidence 
the  anomalous  behavior  of  the  microphone  number  4. 
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Frequency  (Hz) 


Fig.  23.  Power  spectral  density  referred  to  microphone  4  at  Re  =  36  865  computed  for  r= 3.8  (— ),  r= 3.2  ( 


)  and,  r= 2.6  ( . ). 


(p-<p>)/Gp 

Fig.  24.  Probability  density  function  of  pressure  fluctuations  measured  by  microphone  4  at  Re  =  36  865  computed  for  r= 3.8  (□),  r= 3.2  (o)  and  r=  2.6  (A). 

zero  mean  and  unitary  standard  deviation.  The  PDFs  are  all  similar  to  the  reference  Gaussian  curve  except  for  the  case  at 
r= 2.6.  More  specifically,  denoting  with  S  the  Skewness  factor,  it  is  found  S=1.01  for  r  =2.6  and  S=0.06  for  r= 3.8.  The 
origin  of  such  a  behavior  can  be  ascribed  to  the  effect  of  random  pressure  surges,  which  are  statistically  relevant  in  the 
region  of  the  neck,  where  they  might  be  induced  by  the  flow  separation  where  a  significant  pressure  gradient  is  present.  To 
this  extent,  [29]  observed  a  positive  skewness  of  the  wall  pressure  fluctuations  in  the  reattaching  zone  of  a  separated  region 
and  interpreted  this  behavior  as  induced  by  inrush  of  irrotational  flow  towards  the  wall.  A  similar  behavior  has  been  found 
in  the  region  of  flow  separation  generated  upstream  of  a  forward-facing  step  [30].  An  example  of  the  effect  of  these  surges  in 
term  of  wall  pressure  fluctuations  can  be  observed  in  Fig.  25,  where  several  signatures  with  amplitude  over  3ap  can  be 
detected.  To  confirm  that  high  amplitude  pressure  fluctuations  are  related  to  the  recirculation  bubble  dynamic,  the  proper 
orthogonal  decomposition  (POD)  is  applied  to  velocity  field's  domain  T2,  which  is  indicated  in  Fig.  26.  The  POD  is  applied  in 
this  context  to  analyze  the  statistical  properties  of  large-scale  structures  [31],  in  the  present  investigation  generated  within 
the  neck.  When  most  of  the  kinetic  energy  in  the  flow  domain  is  embedded  in  a  large-scale  structure,  a  common  outcome  is 
that  low-order  modes  determine  the  structure  dynamics.  In  this  work  about  1200  snapshots  have  been  used  for  the  POD 
computations  to  ensure  the  convergence  of  the  modes  and  the  achieved  results  will  be  presented  in  the  following.  The  POD 
eigenvalue's  spectrum  is  reported  in  Fig.  27  and  it  is  clear  that  an  important  part  of  the  fluctuating  energy  is  concentrated  in 
the  first  mode.  This  result  seems  to  be  independent  of  r.  As  a  consequence  we  can  calculate  the  first-order  POD  expansion 
taking  as  approximation  of  the  instantaneous  flow  field: 

un  « (un)jv  +  a"Oi 


(10) 


T.  Pagliaroli,  R.  Camussi  /  Journal  of  Sound  and  Vibration  341  (2015)  116-137 


131 


8 

6 

4 


Time  (s) 

Fig.  25.  Time  history  of  pressure  fluctuations  normalized  respect  to  its  standard  deviation  and  measured  by  microphone  4  at  Re  =  36  865  and  r= 2.6. 


Neck 

_ Q 

Cavity 


Fig.  26.  Enlargement  of  the  TVC  geometry  and  POD  domain  representation. 


mode  no. 

Fig.  27.  Spectrum  of  the  eigenvalues  referred  to  r= 2.6  (□)  and  r= 3.8  (A). 


where  the  operator  <  •  )N  represents  the  ensemble  average,  and  the  1st  POD  coefficient  and  mode  respectively. 
Furthermore,  calculating  the  standard  deviation  of  at,  denoted  as  <rfll,  the  following  term  can  be  calculated: 

01=<7ai®,  (11) 

From  statistical  point  of  view,  this  term  can  be  assumed  as  a  model  of  the  most  intense  perturbation  induced  by  the  first 
mode  on  the  mean  field.  Hence,  in  a  probabilistic  sense,  all  instantaneous  velocity  fields  can  be  assumed  to  be  contained  in 
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Fig.  28.  Statistical  reconstruction  of  the  recirculation  bubble  dynamic  using  the  first-order  POD  expansion,  obtained  for  r= 2.6  (a,  b,  c)  and  r=  3.8  (d,  e,  f). 


(a)  (b) 


Fig.  29.  Power  spectral  density  referred  to  microphone  5  obtained  varying  r  and  keeping  Re  constant  (a).  First  mode  frequency,  provided  by  means  of  the 
acoustic  modal  test  (  +  )  and  aeroacoustic  test  (A),  superimposed  to  solutions  of  Eq.  (4)  ( cubical  HR  — )  and  Eq.  (3)  ( double  duct - )  (b). 


the  following  range: 


un  G[(un)N-0i;(un)N  +  0i] 


(12) 


In  Fig.  28,  the  three  fields  corresponding  to  <un)jv-0i,  < un)N  and  < un)N  +  0-i  are  reported  for  the  rectangular  and  the  square 
cavity  (first  and  second  row  respectively).  As  can  be  observed  for  the  geometry  based  on  rectangular  cavity  the  dynamic  of 
the  recirculation  bubble  is  relevant  and  the  large  scale  flow  structure  does  not  seem  to  be  located  within  the  neck.  This 
unsteady  behavior  can  be  the  source  of  the  low  frequency  pressure  surges,  highlighted  previously  by  observing  the  PSDs, 
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(a)  (b) 


2 

Fig.  30.  Power  spectral  density  of  the  wall  pressure  referred  to  microphone  1,  the  pressure  is  normalized  by  qin.  The  data  are  taken  at  Re  =  50  271  and 
r=  2.6,  3.2  and  3.8  (a).  The  data  are  computed  at  Re  =  50  271,  36  865,  30  163  and  r= 3.8  (b). 


Fig.  31.  Normalized  power  spectral  density  referred  to  microphone  1  for  Re  =  50  271  (— ),  Re  =  36  865  ( - )  and  Re  =  30  163  (•  - ).  Spectra  calculated  for 

rectangular  and  square  cavities  are  reported  in  (a)  and  (b).  The  solid  line  represents  the  slope  - 1  whereas  the  dash-dotted  line  the  slope  -2.  Cases  (c)  and 
(d)  report  contour  plots  and  streamlines  of  the  averaged  vorticity  field  in  the  case  of  Re  =  50  271  calculated  for  two  values  of  cavity  aspect  ratio:  r=  2.6  (c) 
and  r=3.8  (d). 
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Fig.  32.  Normalized  power  spectral  density  referred  to  microphone  5  for  different  r  and  Re  =  50  271. 

PDFs  and  the  time  history  of  the  wall  pressure.  As  pointed  out  above,  the  wall  pressure  PSDs  are  influenced  significantly  by 
the  acoustic  mode.  Indeed,  the  acoustic  resonance  affects  the  overall  pressure  dynamics  and  it  is  sensed  by  the  whole  set  of 
microphones  used  in  the  present  experiment.  As  a  matter  of  fact,  the  spectra  reported  in  Fig.  29(a)  clearly  show  a  peak  at 
about  the  same  frequency  as  that  of  the  first  acoustic  mode  determined  by  the  acoustic  numerical  and  experimental 
investigations.  The  acoustic  nature  of  these  peaks  is  demonstrated  by  Fig.  29(b)  where  those  frequencies  are  compared  both 
to  the  analytical  predictions  obtained  by  the  solution  of  Eqs.  (3)  and  (4)  and  experimental  frequencies  measured  by  means  of 
the  modal  test.  According  to  the  acoustic  analysis  of  Section  4,  the  experimental  data  are  in  fairly  good  agreement  with  the 
scaling  law  predicted  by  the  simple  HR  model  and  reported  in  the  following: 


/ 


c0  rr 

2tt  \  rue 


(13) 


The  above  results  clearly  show  that  the  parameter  r  influences  significantly  the  largest  peaks  of  the  pressure  spectra.  Taking 
into  account  this  consideration,  in  order  to  scale  appropriately  the  low  frequency  range  of  the  pressure  spectra,  the 
frequency  can  be  normalized  as  follows: 


strJAmrl/2,  04) 

O' in 

where  Uin  is  the  mean  inlet  velocity  and  din  the  inlet  section  diameter.  Str  is  a  reduced  Strouhal  number  that  takes  into 
account  the  variation  of  r  and,  as  shown  in  Fig.  30,  leads  the  spectra  to  collapse  very  well  in  the  low  frequency  range.  In  the 
plot  of  Fig.  30,  the  y-axis  is  normalized  with  respect  to  the  inlet  dynamic  pressure  (qin)  and  using  the  local  reference  time 
t  =  din/U[n.  The  effect  of  Re  is  reported  in  Fig.  30(b)  where  it  can  be  observed  that  a  collapse  of  the  peaks  is  no  longer 
observed  whereas  the  spectra  collapse  well  in  the  high  frequency  range.  In  order  to  investigate  the  high  frequency  behavior 
at  large  St  (Str  >  0.5)  the  PSDs  for  r=2.6  and  r =3.8,  computed  at  three  different  Re,  are  reported  in  Fig.  31(a)  and  (b).  Two 
different  power  law  decays,  with  exponents  - 1  and  -  2,  are  highlighted  by  the  straight  lines  in  the  log-log  representation. 
The  - 1  power  law  is  observed  at  r=  2.6  and  it  is  characteristic  of  equilibrium  turbulent  boundary  layers  (TBLs),  associated 
with  turbulent  activity  in  the  log  layer  [32  .  On  the  other  hand,  a  power  decay  equal  to  -2  is  typical  of  non- equilibrium  TBL 
and  is  observed  to  apply  in  the  rectangular  cavity  case  (r  =2.6).  Na  and  Moin  [33]  performed  a  DNS  of  a  TBL  developing  over 
a  flat  plate,  and  a  power  law  decay  with  exponent  close  to  -  2  was  observed  for  the  spectra  downstream  of  a  reattachment 
location.  With  reference  to  Fig.  11,  the  APG  at  the  location  of  mic.  1  (that  is  for  x/(L+/)  =  0.46)  is  larger  for  r= 3.8  than  for 
r =2.6.  The  PIV  results  reported  in  Fig.  31(c)  and  (d)  also  support  this  interpretation.  In  fact  it  is  observed  that  the  vorticity 
shear  close  to  the  wall,  in  the  square  cavity  case,  appears  thickened  and  more  irregular  than  in  the  rectangular  one  thus 
suggesting  that  the  wall  flow  in  the  second  case  is  closer  to  a  regular  TBL.  The  pressure  spectra  computed  from  the 
transverse  microphone  (mic.  5)  are  reported  in  Fig.  32  for  different  aspect  ratios.  In  this  case,  the  exponent  of  the  spectral 
scaling  in  the  high  frequency  range  varies  from  -2  to  -4  for  increasing  r.  This  behavior  explained  through  the  vorticity 
fields  provided  by  the  PIV  measurements  is  in  dimensionless  form  coz/core f,  where  core f  =  1/At  and  At  is  the  time  delay 
between  two  consecutive  laser  pulses.  Fig.  31(d)  clearly  shows  that  for  r= 3.8  the  large  scale  vortex  core  is  centered  on  mic. 
5.  Indeed,  as  pointed  out  e.g.  by  [33]  the  power  law  exponent  -4  is  typical  of  a  separation  bubble.  On  the  other  hand,  as 
shown  in  Fig.  31(c),  a  variation  of  the  cavity  aspect  ratio  from  r= 3.8  to  r  =2.6  corresponds  to  a  relevant  displacement  of  the 
vortex  core  thus  explaining  the  power  law  variation  from  -4  to  -2. 
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(a)  (b) 


Fig.  33.  Cross-spectrum  magnitude  ( - )  and  phase  (— )  (a).  Cross-correlation  between  microphones  1  and  2  computed  after  filtering  the  raw  signals 

with  high-pass  filter,  the  curves  are  provided  for  three  different  Re  numbers:  50  271  (— ),  36  865  ( - )  and  30  163  ( —  •)  (b). 


Table  4 

The  convention  velocity  and  the  convection-inlet  velocity  ratio  are  gathered  in  this  table  for  square 
and  rectangular  cavities  at  different  Re  numbers. 


r 

Re 

Uc  (m/s) 

Uc/Uin 

3.8 

30163 

23.6 

0.48 

3.8 

36  865 

29.1 

0.46 

3.8 

50  271 

37.7 

0.49 

2.6 

30163 

24.8 

0.50 

2.6 

36  865 

32.9 

0.52 

2.6 

50  271 

38.0 

0.49 

5.2.  Two-point  pressure  statistics 


The  simultaneous  measurement  of  pressure  by  several  wall  mounted  microphones  allows  for  the  computation  of  cross¬ 
correlations  and  cross-spectra  that  are  crucial  quantities  to  be  analyzed  for  the  characterization  of  the  propagation  of  the 
pressure  perturbations.  Indeed,  the  analysis  of  pointwise  signals  presented  in  the  previous  section  suggests  that  the  low 
frequency  phenomena  are  dominated  by  acoustic  resonances  leading  to  the  generation  of  acoustic  waves  propagating  at  the 
speed  of  sound.  On  the  other  hand,  higher  frequencies  seem  to  be  influenced  by  the  large  and  small  scale  flow  structures, 
localized  in  the  bulk  and  close  to  the  wall  and  inducing  hydrodynamic  pressure  fluctuations.  The  convection  velocity  Uc  is 
directly  dependent  upon  the  nature  of  the  pressure  perturbations.  According  to  e.g.  [34],  it  is  possible  to  define  Uc  through 
the  cross-spectral  phase  <pij  computed  between  a  couple  of  sensors  separated  by  a  distance  £,  as  follows: 


cog 

P'<~Uc{0)g) 


(15) 


The  cross-spectral  magnitude  and  phase,  computed  between  mics.  1  and  2,  is  displayed  in  Fig.  33(a)  for  the  square  cavity 
case.  The  phase  exhibits  a  plateau  at  low  frequencies  and  a  decaying  region  with  a  negative  slope  at  higher  frequencies.  The 
first  one  ranges  from  10  Hz  up  to  the  frequency  of  the  first  acoustic  mode,  the  second  one  ranges  from  the  first  mode 
frequency  up  to  around  2000  Hz.  The  plateau  can  be  ascribed  to  the  acoustic  effect  since  all  pressure  time  scales  propagate 
at  the  same  velocity.  On  the  other  hand,  according  to  the  definition  provided  in  Eq.  (15),  the  negative  slope  can  be 
interpreted  to  be  a  consequence  of  convective  effects.  This  interpretation  is  confirmed  by  the  analysis  of  the  cross¬ 
correlation  between  the  signals  measured  by  mics.  1  and  2  and  properly  filtered  in  advance.  More  specifically,  a  second- 
order  Butterworth  high-pass  filter  using  a  frequency  cut-off  of  200  Hz  has  been  applied.  All  signals  are  processed 
synchronously,  to  avoid  the  occurrence  of  phase  drift  across  the  records  from  different  microphones.  In  this  way,  the 
acoustic  resonance  is  excluded  and  only  high  frequencies,  that  are  supposed  to  be  associated  to  the  hydrodynamic  pressure, 
are  considered.  The  cross-correlations  between  the  filtered  signals  are  reported  in  Fig.  33(b),  where  can  be  observed  that  the 
time  delay,  t,  associated  to  the  peak  of  the  correlations  vary  with  Re.  This  result  according  to  the  cross-spectra  phase  can  be 
ascribed  to  the  hydrodynamic  advection.  This  analysis  allows  us  the  computation  of  the  hydrodynamic  convection  velocity, 
Uc  =  £/r,  whose  values,  for  different  r ,  are  summarized  in  Table  4.  It  can  be  observed  that  the  ratio  between  Uc  and  the 
velocity  inlet  is  about  constant  and  equal  to  0.5,  furthermore  it  is  not  affected  neither  by  r  nor  by  Re. 
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6.  Concluding  remarks 

The  acoustic  and  the  aeroacoustic  behavior  of  a  RPE  have  been  studied  by  applying  several  experimental  measurement 
techniques  and  an  acoustic  numerical  simulation.  The  fluid  dynamic  field  for  different  geometrical  configurations  has  been 
characterized  using  the  PIV  technique.  The  aerodynamic  study  evidenced  the  main  flow  structures,  including  separations  bubbles 
near  the  walls  and  a  large  scale  vortex  in  the  bulk  and  it  has  been  verified  that  the  square  cavity  geometry  is  the  best  in  terms  of 
large  scale  vortex  intensity.  From  the  acoustic  view  point  it  has  been  highlighted  that  RPE  behaves  as  a  Helmholtz  resonator,  and  a 
novel  end-correction  formula  has  been  proposed  and  validated  against  experimental  data.  Further  a  prediction  model  based  on 
Helmholtz  Resonator  formulation  and  on  a  novel  end-correction  has  been  derived.  The  aeroacoustic  measurements  have  shown 
that  the  first  acoustic  mode  dominates  the  wall  pressure  spectra  also  in  the  presence  of  a  mean  flow.  A  reduced  form  of  Strouhal 
number  has  been  proposed  on  the  basis  of  the  frequency  scaling  of  the  first  acoustic  mode  and  a  good  collapse  of  the  frequency 
spectra  in  low  frequency  range  has  been  achieved.  On  the  other  hand,  it  has  been  shown  that  the  high  frequency  region  of  the 
spectra  is  affected  by  the  APG  that  influences  the  jet-wall  interaction  and  induces  a  non-equilibrium  condition  of  the  flow  close  to 
the  wall.  This  behavior  induces  a  variation  of  the  high  frequency  spectra  slope,  from  - 1  to  -  2.  The  recirculation  region  has  been 
characterized  also  through  a  POD  decomposition  of  the  PIV  velocity  fields.  It  is  found  that  the  dynamics  of  the  bubble  formed  into 
the  neck  depends  on  r,  leading  the  pressure  fluctuations  to  increase  at  low  frequencies.  The  results  provided  in  the  present 
investigation  definitely  demonstrate  that  the  main  parameter  to  be  taken  into  account  in  the  design  of  this  kind  of  geometry  is  the 
cavity-neck  aspect  ratio. 
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